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Abstract 



We prove that solutions of a mildly regularized Perona- Malik equation converge, in 
a slow time scale, to solutions of the total variation flow. The convergence result is 
global-in-time, and holds true in any space dimension. 

The proof is based on the general principle that "the limit of gradient-flows is the 
gradient-flow of the limit" . To this end, we exploit a general result relating the Gamma- 
limit of a sequence of functionals to the limit of the corresponding maximal slope curves. 
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1 Introduction 



The Perona-Malik equation 



u t = div 




(1.1) 



(where 11 C M" is a bounded open set) is arguably the most celebrated example of 
forward-backward diffusion process. It was introduced by P. Perona and J. Malik [34] 
in the context of image denoising. It is formally the gradient-flow of the functional 



The forward-backward nature of (1.1) depends on the convex- concave behavior of 
the integrand in (1.2). Equation (1.1) has generated a considerable literature (see [3, 
7, 8, 9, 10, 11, 13, 14, 15, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 37, 38]), focussed both 
on numerical and on analytical aspects. The big challenge is to reconcile the empirical 
and practical efficacy of the method, supported by several numerical observations, with 
the expected analytical ill-posedness of a backward parabolic equation. A satisfactory 
theory would represent a solution to the Perona-Malik paradox, as named after [31], but 
it still seems to be out of reach. 

In a recent paper, P. Guidotti [30] introduced a mild regularization of (1.1). He 
considered the family of functionals 



where 5 > is a parameter. The integrand remains nonconvex, at least in the interesting 
cases where S is small, but now it is convex-concave-convex, and it grows quadratically 
at infinity. This is enough to guarantee that the corresponding gradient-flow equation, 
which simply reads as 



has a unique global-in-time solution in the sense of Young measures, according to the 
theory developed in [32, 19]. 

Several qualitative properties of these solutions are reported in [30]. In particular, 
it seems that the well-known staircasing effect of the original Perona-Malik equation is 
now replaced by a "ramping" effect, namely the tendency of solutions to alternate flat 
plateaus and bounded growth regions in a piecewise fashion (see all figures in [30]). 

In [30] it is also observed that the limit of approximated solutions as 8 — > + is 
the trivial stationary solution frozen in the initial condition, and "thus the only way to 
produce a meaningful limit would involve a rescaling of time in the process" . 




(1.2) 





(1.3) 
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In this paper we follow this path. We take a family of solutions vs(t) of the regularized 
model (with initial conditions vs(0) converging to some u ), and then we speed up the 
evolution by considering the family of functions 

u £ (t) := v 4 -i £ 2 llog£l t > 0- (1-4) 

This happens to be the right rescaling factor, in the sense that u £ (t) uniformly 
converges to a nontrivial limit u(t) as e — > + . Moreover u(t) turns out to be the 
solution of the total variation flow (see [6, 4]) 

u t = div (i,t)efix(0,+oo), (1.5) 

with initial condition uq. Any slower time rescaling leads in the limit to the stationary 
solution u(t) = Mo, and any faster time rescaling produces a limit which for every t > 
coincides with the constant function equal to the average of uq in Q. 

All qualitative properties of the total variation flow (see [5]) are consistent with the 
numerical experiments presented in [30]. Therefore, our convergence result provides a 
rigorous justification of all these properties. 

The proof of the convergence result involves three main steps. First of all, we inter- 
pret u £ (t) as a gradient-flow. As it comes from [30], u e (t) is the solution in the sense of 
Young measures of the forward-backward equation associated to the formal gradient-flow 
of the nonconvex energy 

E £ (u) := / <p £ (\Vu(x)\) dx, (1.6) 
Jn 

where 

The key point is that u e (t) is also the gradient-flow of the relaxed (convex) energy 

E;*(u):= [ ^r(|V«(x)|) dx, (1.7) 
Jn 

where f** is the convexification of ip £ . 

In other words, as a result of the first step we can forget about Young measures and 
forward-backward equations, and think of u £ (t) as the solution of a degenerate forward 
parabolic equation, or better as the gradient-flow of a convex (although not strictly 
convex) functional. We stress that this is a general fact. Solutions provided by the 
theory developed in [32, 19] always coincide with gradient-flows of the corresponding 
convexified (or relaxed) energies. Quite surprisingly, this has never been observed in the 
literature up to our knowledge. 
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In the second step we compute the Gamma-limit of the energies (Theorem 2.1), and 
we discover that 



where TV(u) denotes the total variation of u. Since the total variation flow is the 
gradient-flow of TV(u), our convergence result is now equivalent to say that the limit 
of gradient-flows is the gradient-flow of the Gamma-limit. 

This is the content of the third step, where we deduce it from a general result 
(Theorem 3.1) which we state and prove in the abstract setting of maximal slope curves 
in metric spaces (see Section 3 for further details and references). We believe that the 
scope of the general result goes far beyond this simple application. 

We conclude by discussing how this problem relates to recent investigations about 
the "slow time" behavior of approximations of the Perona-Malik equation. The starting 
point is the observation that the evolution, despite of the different approximation meth- 
ods, seems to exhibit always three different time scales, named "fast time", "standard 
time" , and "slow time" according to [8] . 

In a fast time of order o(l) solutions develop microstructures in the concave regime 
(staircasing). In a time scale of order 0(1) (standard time) solutions behave as expected 
in the original model, with a smoothing effect in the concave regime, and sharpening 
of regions where the gradient is large. At the "end" of standard time, solutions have 
an almost piecewise constant structure, and this is consistent with the intuitive idea 
that piecewise constant functions are stationary points of PM(u). On the other hand, 
only constant functions (and not piecewise constant functions) are stationary points of 
the usual approximating models. As a consequence, approximating solutions exhibit a 
transition from a piecewise constant structure to a constant value (equal to the average 
of Uq). The transition turns out to be very slow because there is almost no energy left. 

This gives rise to the so called "slow time" motion, in which the plateaus of the 
piecewise constant function move in the vertical direction, with jump points which 
remain fixed in space. Although the existence of this phase seems to be independent of 
the approximation method, the law of the vertical motion does depend on it. 

All previous results on this problem are limited to the one dimensional case. Let 
it be a piecewise constant function defined in an interval, and let S u be the (finite or 
countable) set of its jump points. Let J x denote the jump height in a point x G S u , and 
let us consider the following energies 



In the case of a fourth order regularization of (1.1), corresponding to adding a vanish- 
ing second order term to (1.2), G. Bellettini and A. Fusco [8] conjectured that the vertical 
motion is governed by the gradient-flow of H 1 / 2 (u). They supported their conjecture by 
proving the corresponding Gamma-limit result for the energies. The missing step is a 





H a (u) := W (with a E (0, 1]), H (u) := ^ log \J X 
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rigorous proof that also in that case the limit of gradient-flows is the gradient-flow of 
the limit. Unfortunately Theorem 3.1 does not apply to their functionals. 

In the case of the semidiscrete scheme (see [10, 15]) the vertical motion is governed 
by the gradient-flow of H (u), which in a certain sense represents a limit case. The proof 
given in [15] exploits a variant of Theorem 3.1, complicated by the fact that H (u) is 
not bounded from below. 

What we show in this paper is that in the model proposed in [30] the slow-time 
vertical motion is governed by the gradient-flow of Hi{u). This is the opposite limit 
case, and phenomena are completely different. The good news is that the limit energy 
is convex. This simplifies the analysis, which here can be carried out in any space 
dimension, and delivers a well known limit problem. On the other hand, the relaxation 
of Hi{u) is TV(u), hence it is finite in the whole space of bounded variation functions. 
As a consequence, in this case the slow time motion is not limited to piecewise constant 
functions. This is hardly surprising after reminding that in this model the motion in 
standard time is trivial. 

This paper is organized as follows. In Section 2 we fix notations and we state our 
convergence results. Section 3 is devoted to limits of maximal slope curves in metric 
spaces. In Section 4 we prove our main result. 

2 Notations and statements 

Let n be a positive integer, and let Q C R™ be an open set, which for simplicity we 
assume to be bounded and an extension domain (see Definition 3.20 in [1], satisfied 
by all bounded open sets with Lipschitz boundary). The more general ambient space 
we consider is L 2 (Q). We write or simply \\u\\ p , to denote the p-norm (with 

p G [l,+oo]) of a function u G L P (Q). All the energies we consider, and in particular 
E £ {u) and E**(u), are always thought as defined in the whole space L 2 (fl) by setting 
them equal to +oo outside their natural domain. We write BV(Q) to denote the space 
of all functions u G L 2 (Q) with finite total variation TV(u). Once again, we think of 
TV(u) as defined for every u G L 2 (Q), with TV(u) < +oo if and only if u G BV(Q). 

For every 5 > we consider equation (1.3), with Neumann boundary conditions, 
and an initial datum. It has been shown in [30] that solutions v$(t) exist in a suitable 
weak sense. For every e G (0, 1), we define u £ (t) by rescaling v$(t) according to (1.4). It 
turns out that u £ (t) is a solution in the same weak sense of equation 





(2.1) 



with Neumann boundary conditions 



du £ 



(x,t) = 



(x,t) G <9f2 x (0, +oo), 



(2.2) 



On 



4 



and initial condition 

u e (x, 0) = Uo £ (x) x G Q. (2.3) 
In the following result we collect properties of u £ {t). 

Theorem A (Properties of rescaled approximating solutions) Let n be a posi- 
tive integer, let Q C R™ be a bounded extension domain, let e G (0, 1), and let u 0s G 
L 2 (f2). Let E** be the functional defined in (1-7). 
Then the following properties hold true. 

(1) (Weak Young measure solution and regularity) There exists a unique function u £ {t) 
and a (not necessarily unique) gradient Young measure u £ inQx [0, +oo) such that 
the pair (u £ , u £ ) is a weak Young measure valued solution of problem (2.1) through 
(2.3) in the sense of [32, 19]. 

Moreover, we have that 

u £ G C° ([0, +oo); L 2 (Q)) n C 1 ((0, +oo); L 2 (Q)) , 

and for every t > we have that u £ (t) is regular enough so that the right-hand 
side of (2.4) lies in L 2 (Q). 

(2) (Degenerate forward parabolic equation) The function u £ of statement (1) is the 
unique solution in Q x (0, +00) of the partial differential equation 



u £t = -VE* £ *(u £ ) =div 



(vC)'(lv^l) V "- 



|Vu e | 

with Neumann boundary conditions (2.2), and initial condition (2.3). 



(2.4) 



(3) (Gradient-flow integral inequality) The function u £ of statement (1) is the unique 
function satisfying (2.3) and the inequality 

E:*(u e (s))-E* £ *(u e (t)) > ^'llu^llldr+i jf'llV^'Mr^llidr (2.5) 
for every < s < t. 

(4) [L v estimate) If u 0e G L P (J7) for some p G [l,+oo], then u s (t) G L p (Vl) for every 
t > 0, and the function t — > \\u £ (t)\\LP(n) is nonincreasing. 

Theorem A above shows that the function u e (t) can be characterized in at least 
three different ways, either as the solution in the sense of Young measures of a forward- 
backward equation, or as the solution of a forward degenerate parabolic equation (for 
example in the sense of [12]), or as a maximal slope curve (gradient-flow inequalities). 
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What we need in this paper is only the last one. Solutions generate a contraction 
semigroup in L 2 (fl). 

As far as we know this equivalence has never been stated explicitly. Nevertheless, 
it follows from some general facts which nowadays are quite well known, and which we 
now recall briefly. 

First of all, the three approaches lead to a unique solution. Uniqueness follows in the 
first case from the strong requirements imposed on the structure of the corresponding 
Young measure, in the second case from the contraction property of the semigroup 
generated by a forward parabolic equation, in the third case from the convexity of the 
energy E**(u). 

Secondly, in all three approaches the solution is usually obtained (or at least it 
can be obtained) as the limit of approximated solutions constructed via an iterated 
minimization process, known as minimizing movement (see [17]). As already observed 
in [32, 19, 30], the minimization procedure gives the same result when applied to E £ {u) 
or E**(u). 

In conclusion, all three approaches define a unique solution through an analogous 
procedure, hence the solution is the same. 

In this paper we are interested in the behavior of u £ (t) as e — > + . Following the 
gradient-flow approach, the first thing to do is understanding the limit behavior of the 
energies. This is the content of next result. 

Theorem 2.1 (Gamma-convergence and compactness) Let n be a positive inte- 
ger, and let Q C ]R n be a bounded extension domain. For every e > 0, let E e : L 2 (Q) — > 
[0, +oo] and E** : L 2 (Q) — > [0, +oo] be defined as in (1.6) and (1.7), respectively, if 
u G H 1 ^), and +oo otherwise. 

Then we have the following conclusions. 

(1) (Gamma-convergence) We have that (1.8) holds true with respect to the topology 

fL 2 (n). 

(2) (Compactness) Let {w e } eG (o,i) ^ L 2 (Q) be a family of functions such that 

sup {KIU + E**(u e )} <+oo. (2.6) 

ee(0,l) 

Then the family {u £ } is relatively compact in L 2 {Vl). 

The gradient-flow of the limit functional TV{u) is the so called total variation flow, 
and the corresponding partial differential equation is (1.5). The right-hand side of (1.5) 
needs to be interpreted in a suitable weak sense when the gradient vanishes, and this 
happens in large regions because solutions tend to develop flat plateaus. A general 
existence and uniqueness result was proved by F. Andreu, C. Ballester, V. Caselles, 
and J. M. Mazon [4] (see also [5]) using the theory of accretive operators in Banach 
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spaces. In [4] the operator in the right-hand side of (1.5) is interpreted as the limit of 
the p-Laplacian as p — > 1 + . 

In this paper we limit ourselves to initial data in L 2 (Q), in which case existence of a 
unique solution is provided also by the theory of maximal monotone operators [12], as 
explained in [6]. In this context the right-hand side of (1.5) is the sub differential of the 
convex functional TV(u). As in the case of approximating problems, this formulation 
is equivalent to the gradient-flow integral inequality 

TV(u(s)) - TV(u(t)) > \ jf |K(r) Wldr + lJ* || VTV(«(r))||| dr 

for every < s < t, where || VTV(u(t))\\ 2 is the minimal norm of an element in the 
sub differential of the functional TV in the point u(r), which in turn coincides with the 
slope of TV in u{r) as defined in Section 3 in an abstract metric setting. This is the 
characterization of the total variation flow which we need in this paper. 

Our main result is the convergence of u E (t) to the solution of the total variation flow 
with the same boundary conditions, and the limit initial datum. We point out that we 
do not assume initial data to be a recovery sequence, and we do not ask their energy to 
be bounded. 

Theorem 2.2 (Global-in-time convergence) Letn be a positive integer, letVt C W 1 
be a bounded extension domain, let u G L 2 (Q), and let {ttoe}ee(o,i) — L 2 (Q) be a family 
of functions such that 

lim uq £ = uo in L 2 {Vt). (2.7) 

For every e G (0, 1), let u e be the solution of the rescaled approximating problem with 
initial condition uq £ , in the sense of Theorem A. Let u{t) be the solution of the total 
variation flow with Neumann boundary conditions and initial datum uq. 

Then we have that u £ {t) — > u(t) in C° ([0, +oo); L 2 (Q)), namely 

lim sup \\u £ (t) - u(t)\\ L 2(n) = 0. (2.8) 

We conclude this section with a heuristic argument which justifies the rescaling 
leading to (1.6). 

Let us consider the one dimensional case where Q = (—1,1), let J > 0, and let 
v E BV({—\, 1)) be the piecewise constant function equal to — J/2 in (—1, 0), and equal 
to J/2 in (0, 1). Let rj > 0, and let h G if^R) be a function such that h(x) = J/2 for 
every x > rj, h(x) = — J/2 for every x < —i], and h'(x) > for almost every x G (-T), rj). 
For every e > 0, let v £ (x) be the approximation of v defined as v £ (x) := h(x/e) for every 
■ri i l.D. 

For every e < r/^ 1 , plugging v £ into (1.6), with a variable change we obtain that 
E £ (v £ ) = £ log ^1 + 1 [h'(x)] 2 ^ dx + \ £ [h\x)f dx. 
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Letting e — > + we find that 

1 P 9 

liminf E £ {v £ ) > 2r] + - / [h'{x)\ dx. (2.9) 

In order to estimate the right-hand side, we first minimize with respect to h, and 
we discover that the optimal choice is the function h(x) defined as J(2r]) _1 a; for every 
x G (—77, rj). Then we compute the integral, and finally we apply the inequality between 
arithmetic and geometric mean to deduce that 




(2.10) 



Estimates (2.9) and (2.10) suggest that the cost of a jump could be the jump height, 
which leads to conjecture that the Gamma-limit is TV{u) for a general u in any dimen- 
sion. 



3 Passing to the limit in maximal slope curves 

The abstract theory of gradient-flows in metric spaces was introduced in [16], and then 
developed by the same authors and collaborators in a series of papers (see [18, 33] and 
the references quoted therein). For a modern presentation we refer to [2]. Here we just 
recall some basic definitions. 

Let (X, d) be a metric space, let R := MU {— 00, +00} be the extended real line, and 
let F : X — > R be any function. The (descending) slope \VF\(x) of F in x is defined to 
be +00 if F(x) G" K, and otherwise 

lr7nl , x .. max\F(x) - F(y), 0} r ^ . 

VF \(x) := hmsup \ e P, +00]. 

y^x d{-'-. !l) 

For every T > 0, the space AC 2 ([0, T]; X) is the set of all functions v : [0, T] ->■ X 
for which there exists g G L 2 ((0,T)) such that 

d(v(t),v(s)) < J g(r)dT V0 < s < t < T. (3.1) 

It can be seen that there exists a smallest function g(t) satisfying (3.1). This function 
is called the metric derivative of v, and it is denoted by |f'|(t). 
A maximal slope curve for F in [0,T] is a triple (u,ip,E) where 

• tie AC 2 ([0,T};X), 
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• ip : [0,7"] — > R is a nonincreasing function such that for every < s < t < T we 
have that 

m - m >±J* \uf(r) dr+^J* |VF| 2 («(r)) dr, (3.2) 

• EC [0,T] is a set with Lebesgue measure equal to such that 

^(f) = F(u(f)) Vt G [0, T] \ £. (3.3) 

To be more precise, the second integral in the right-hand side of (3.2) should be 
an upper integral, since at this level of generality there is no reason for the function 
t — > \VF\(u(t)) to be measurable. On the other hand, it can be easily proved that 
actually it is always true that \u'\(t) = \VF\(u(t)) for almost every t G [0,T], which 
implies the required measurability. 

When F is a C 1 function in a Hilbert space X, this weak formulation is equivalent 
to the classical one, namely to asking that u'(t) = —VF(u(t)) for every t G [0, T\. 

Besides generality, the advantage of this weak formulation is that inequalities and 
integrals are more stable than equalities and derivatives. It follows that maximal slope 
curves exist under general assumptions on F (see Theorem 2.3.1 in [2]), and are quite 
stable when passing to the limit, both with respect to initial conditions, and with respect 
to functionals. Results in this direction are contained in [18] and [35] in a Hilbert setting, 
and in [36] in a metric setting, but assuming that initial data are a recovery sequence. 

Here we state a quite general result, used in a special case also in [15]. 

Theorem 3.1 (Limits of maximal slope curves) Let X be a metric space, let F : 
X — > R be a function, and let F n : X — >■ R be a sequence of functions. 

Let us assume that for every x G X, and every pair of sequences {n^} C N and 
{xk} CI, we have the implication 



Uk — > +00, Xk — > x 

sup {\F nk (x k ) \ + \ VF nk \(x k )\ < +oo 
ken 



Jjm^F nk (x k ) = F(x). 
liminf \VF nk \{x k ) > \VF\{x). 

k— >+oo 



(3.4) 



Let T > 0, and for every n G N let (u n , ip n , E n ) be a maximal slope curve for F n in 
[0, T] . Let us assume that 

sup max \ib n (t)\ < +oo, (3.5) 

n&i te[0,T\ 

and that {u n } has a pointwise limit, namely there exists u : [0,T] — > X such that 

lim u n (t) = u(t) Vt G [0,T]. 

Then there exist ip and E such that (u, ip, E) is a maximal slope curve for F . 



9 



Proof From the definition of maximal slope curve, for every n G N we have that 

Ms) - Mt) >\f K\ 2 (r) dT + ^£ \VF n \ 2 (u n (r)) dr (3.6) 

for every < s < t < T, and 

Mt) = F n (u n (t)) V* G [0, T] \ E n . (3.7) 

The functions ip n {t) are nonincreasing, and equi-bounded because of (3.5). Thus 
the usual compactness result for monotone functions (known as Helly's Lemma, see for 
example [2, Lemma 3.3.3]) implies the existence of a nonincreasing function : [0, T] — > 
[0, +oo) such that (up to subsequences, not relabeled) 

lim = Vt€[0,T]. (3.8) 

n— >+oo 

Let us consider now the right-hand side of (3.6). Setting s = and t = T, and using 
once more assumption (3.5), we obtain that 

sup / \u' n \ 2 (t) dr < +oo, (3.9) 
nen Jo 

sup/ \V F n \ 2 (u n (r)) dr < +oo. (3.10) 

neN Jo 

From (3.9) we easily deduce that u G AC 2 ([0,T];X), and 

liminf / |<| 2 (r)rfr> / |M| 2 (r)rfr V0 < s < t < T. (3.11) 
From (3.10) and Fatou's Lemma we obtain that 

/ ( liminf |VF„| 2 (u n (r)) \ dr < liminf / |VF n | 2 (u n (r)) dr < +oo, 

Jo \ n ^+°° J n^+oo J Q 

hence there exists a set E' C [0, T], with Lebesgue measure equal to 0, such that 

liminf |VF n | (u n (t)) < +oo Vt G [0,T] \ E'. (3.12) 

n— »+oo 



Let us introduce the set 



E:=E'U I |J , 

VneN / 



which has clearly Lebesgue measure equal to 0. Let us consider any t G [0, T] \ E. Since 
t ^ from (3.7) and (3.5) we have that 

sup \F n (u n (t))\ = sup \ip n (t)\ < +°°- ( 3 - 13 ) 

neN neN 
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Moreover, due to (3.12) there exists a (t-dependent) sequence n k — > +00 such that 
lim \VF nk \(u nk (t)) = liminf \VF n \(u n (t)) < +00. (3.14) 

k— >+oo n— »+oo 

Thanks to (3.13) and (3.14), sequences {n k } and {u nk (t)} satisfy the assumptions in 
the left-hand side of (3.4). If follows that for every t G [0,T]\E we have that 

V(t)= lim Mt)= lim Vn fc (0 = lim ^» fc K fc (*)) = 

n— >+oo fe— >+oo fc— >+oo 

which proves (3.3), and 

|VF|(u(t)) < liminf |VF n JK fe (t)) = lim | VF nk \(u nk (t)) = liminf \V F n \(u n (t)) , 

fc->+oo fe^+oo 

so that one more application of Fatou's Lemma gives that 

/ |VF| 2 {u{r))dr< liminf / | VF„| 2 (u„(r)) rfr VO < s < t < T. (3.15) 

We can now take the liminf of both sides of (3.6). Thanks to (3.8), (3.11), and (3.15) 
we obtain that (3.2) holds true. This completes the proof that {u,ip,E) is a maximal 
slope curve for F. □ 

Thanks to Theorem 3.1 above, any convergence result for maximal slope curves is 
reduced to verifying three assumptions. The first one is the existence of a pointwise 
limit, namely a compactness result. The second one is estimate (3.5), which in general 
follows from suitable assumptions on the sequence of initial data and some boundedness 
from below of the functionals. The third and more important assumption is (3.4). In 
the last part of this section we show that (3.4) follows from Gamma-convergence in a 
class of functionals which contains all convex functionals in Banach spaces. 

Definition 3.2 (Slope Cone Property) Let (X, d) be a metric space. A function 
F : X — > M satisfies the Slope Cone Property if 

F(y)>F{x)-\VF\(x)-d(x,y) 

for every y G X and every x G X such that F(x) G R and |VF|(x) < +00. 

Remark 3.3 Let X be a Banach space. Then every convex function F : X — > [0, +00] 
fulfils the Slope Cone Property. If in addition F is lower semicontinuous, then its slope 
|VF|(x) is lower semicontinuous. 

On the other hand, also in Banach spaces there do exist nonconvex functions satis- 
fying the Slope Cone Property. An example is F(x) = —\x\ in R. 

Proposition 3.4 Let X be a metric space, and let F n : X — >■ R be a sequence of 
functions satisfying the Slope Cone Property. Let us assume that there exists 

F(x) :=T- lim F n (x). (3.16) 
Then the sequence {F n } satisfies assumption (3.4) of Theorem 3.1. 
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Proof Let n k — > +00 and x k — > x be two sequences as in assumption (3.4). Let M be 
the supremum in the left-hand side of (3.4). Let z n — > x be a recovery sequence for x, 
namely a sequence such that F n (z n ) — > F(x). 

From the Slope Cone Property, and the uniform bound on slopes, it follows that 

Fn k {z nk ) > F nk (x k ) - \V F nk \(x k ) ■ d(x k , z nk ) > F nk (x k ) - M d(x k ,z nk ). 

Taking the lim sup of both sides we obtain that 

F(x) = lim sup F nk (z n J > lim sup F nk (x k ). 

The opposite inequality with the liminf follows from assumption (3.16). This proves 
the first limit in the right-hand side of (3.4). 

Let us prove now the liminf inequality for slopes. Let L denote the liminf in the 
left-hand side of (3.4). Let us take any y £ X, and let y n — > y be a corresponding 
recovery sequence. Due to the Slope Cone Property we have that 

Fn k {Un k ) > F Hk 

(xk) - \VF nk \(x k ) ■ d(x k 

) Vn k ) ■ 

We already proved that the first term in the right-hand side tends to F(x). Therefore, 
taking the lim sup of both sides we obtain that F(y) > F(x) — Ld(x,y). Since y is 
arbitrary, this easily implies that L > \VF\(x), which completes the proof. □ 

We conclude by mentioning a straightforward extension of Definition 3.2 and Propo- 
sition 3.4 in the same spirit of [18, 33]. 

Remark 3.5 One can weaken Definition 3.2 by asking that F : X — > 1U {+00} (so we 
exclude the value —00), and there exists a continuous function $ : X 2 xR 3 4l such 
that x, u, v, w) — for every (x, u, v, w) £ X x IR 3 , and 

F(y) > F(x) - |VF|(x) • d(x,y) - $ (x,y, F(x), F(y), \VF\(x)) ■ d(x,y) 

for every (x,y) £ X 2 such that F(x) £ R, F(y) £ R, and \VF\(x) £ R. 

We call this property "<3> Slope Cone Property" . In a Banach space it is fulfilled, for 
example, by the sum of a C 1 function and a convex function. 

It can be easily proved that Proposition 3.4 holds true also if all functions F n satisfy 
the $ Slope Cone Property with respect to the same function $. 

4 Proofs 

4.1 Gamma-convergence and compactness 

The first equality in (1.8) is a general property of Gamma-convergence. So we can 
concentrate on the second one. 
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A standard approach, suggested by the heuristic argument at the end of Section 2, in- 
volves a blow up argument in order to reduce the Gamma-liminf inequality to minimizing 
the right-hand side of (2.9) with respect to r\ and h, and a proof of the Gamma-limsup 
inequality via the density of piecewise constant functions with smooth level sets, for 
which a recovery sequence can be constructed by adapting the optimal profile in the 
direction orthogonal to level sets. 

In both cases we follow a different and more elementary approach, which exploits 
that our functionals depend only on the gradient. 



Gamma-liminf inequality We claim that for every a G (0, 1) and b G (0, 1) there exists 
E\ G (0, 1) such that 

<p**(a) > a\a\-b \/a G R, Ve G (0, £\). (4.1) 
If we prove this claim, then for every u G L 2 (Vt) and every e G (0, £i) we have that 

E**(u)= [ ip** (\Vu(x)\) dx > [ (a\Vu{x)\-b)dx = aTV{u)-b\Vt\. 
Jn Jn 

Since the functional TV{u) is lower semi continuous, this proves that 



r-liminf E**(u) > aTV{u) - b\Q\ 



Letting a — > 1~ and b — > + , we obtain the required inequality. 

So we are left to prove (4.1). Since the right-hand side is convex, it is enough to 
prove (4.1) with (p £ instead of </?**. Moreover, without loss of generality we can assume 
that a > 0. Now we distinguish four cases. 

If a G [0, b], then we have that 

{ Ps(< J ) > > ab — b > aa — b. 
If a G [b, (e 2 - 1) 1/2 ], then we have that 

1 



<p £ (a) -aa> — log(l + b 2 ) - aVe 2 - 1. 

2e\ ioge\ 

The right-hand side tends to +oo as e — > + , so it is greater than —b when e is small 
enough. 

If a G [(e 2 - 1) 1/2 , (e\ \oge\y 1 ], then we have that 

1 2 a I- a 

(fey®) — a<J > —— 1 loge 



2s | log ^ | e|loge| e|loge|' 

so that the conclusion follows as in the previous case. 
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Finally, when a > (e\ loge|) 1 we apply the inequality between arithmetic mean and 
geometric mean, and we obtain that 

log a e 2 ( log a \ 1/2 / 1 , ( 1 x ^ "'" 
<Pe{?) > -r, r + t 0" > c • n r > c • i n r log 



e|loge| 4 \ | log e | y '" [|log£| \e|loge 

It is not difficult to see that the coefficient of a tends to 1 as e — > + , hence it is 
greater than a when e is small enough (in this point it is essential that a < 1). 
This completes the proof of (4.1), hence also of the Gamma-liminf inequality. 



Gamma-limsup inequality By a classical density argument it is enough to find a recov- 
ery sequence for all functions u G C l {Vl) whose gradient is bounded in Q. To this end, 
it is enough to show that for any such function we have that 

limsup / <p** (\Vu(x)\) dx < TV{u). (4.2) 

In turn, (4.2) is proved if we show that for every M > we have that 

<p?(<r)<a e \<r\ Va G [-M, M] 

for a suitable coefficient a £ which tends to 1 as e — > + . Let us assume, without loss 
of generality, that a G [0, M]. Since ip** is the convexification of ip E , we can estimate 
(fil*(cr) from above with the linear function interpolating the values of (p £ in and 2e _1 . 
As soon as M > 2s~ 1 we obtain that 

As required, the coefficient of a in the right-hand side tends to 1. This completes 
the proof of (4.2). 

Compactness From (4.1) with a = b = 1/2 we obtain that 

ET{Ue) -li (|VMe(;r)l ~ 1 ) dx = \ TV ^) 

This estimate and assumption (2.6) yield a uniform bound on the L°°-norm and on 
the total variation of u £ . Thanks to well known embedding theorems in BV(Q) (see [1], 
this is the point where we need Q to be an extension domain), this implies that the 
family {u £ } is relatively compact in L P (Q) for every p < +oo, and in particular in 
L 2 (n). □ 
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4.2 Convergence of approximating solutions 

In the first three paragraphs we prove the result with the further assumption that 

sup (IMIco + ^rKte)} < +°°- ( 4 - 3 ) 

e€(0,l) 

Then in the last paragraph we prove it for general data. 

Compactness on bounded time intervals We show that, for every T > 0, the family 
{u £ } is relatively compact in C° ([0, T]; L 2 (Q)). This follows from Ascoli's theorem 
provided that we show that solutions are l/2-H61der continuous with equi-bounded 
Holder constants, and that for every fixed t > the family {u £ (t)} C L 2 (Vt) is relatively 
compact. 

From (2.5) and Holder's inequality we have that 

\\u £ (t) - u e (s)h < f s H(r)\\ 2 dr < \t-s\^y\\u' £ (r)\\ldr^ ' 

< \t-s\ 1 / 2 {2E;*(u 0£ )} 1/2 , 

so that the uniform bound on Holder constants follows from assumption (4.3). 

Moreover, from (2.5) and statement (4) of Theorem A we have also that the functions 
t — > E**(u £ (t)) and t — >■ ||u e (i)||oo are nonincreasing, hence 

IMOIloo + E**(u £ (t)) < lluoelU + E;*(u 0£ ) Vt > 0. (4.4) 

Thanks to assumption (4.3), the right-hand side is bounded independently of e. 
Therefore the compactness result in statement (2) of Theorem 2.1 implies that the 
family {u £ (t)} is relatively compact in L 2 (f2) for every fixed t > 0. 

Characterization of the limit We prove that, for every T > 0, any limit point in the 
interval [0,T] of the family {u £ } of approximating solutions is the solution u(t) of the 
total variation flow in the same interval. Since the solution of the limit problem is 
unique, this is enough to prove the convergence of the whole family. 

Let e n — > + be any sequence such that u £n (t) uniformly converges to some v (t) in 
[0,T]. By (2.7) we have that v (0) = u . So it is enough to show that v (t) is a maximal 
slope curve for the functional TV(u) in [0, T\. 

To this end, we apply Theorem 3.1 to the sequence of functionals {E**(u)}. Indeed 
they are convex functionals, and they Gamma-converge to TV(u) because of Theo- 
rem 2.1. Thus Remark 3.3 and Proposition 3.4 prove that assumption (3.4) of The- 
orem 3.1 is satisfied. Also (3.5) holds true because of (4.4), and the fact that the 
functionals are nonnegative. 

Since all the assumptions are satisfied, Theorem 3.1 implies that v (t) is a maximal 
slope curve for the functional TV(u). 
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Uniform convergence for all positive times It remains to prove that the convergence is 
global-in-time, as stated in (2.8). This follows from two general facts. 

The first one is that u(t) tends, as t — > +00, to the constant function m m equal to 
the average of u in Q (see [5]). The second fact is that u £ (t) — is the solution of the 
approximating problem with initial datum uq £ — Moo, hence statement (4) of Theorem A 
(with p = 2) implies that t — > ||ii £ (i) — W00II2 is a nonincreasing function. 

Therefore for every T > we have that 

SUp \\u £ (t) - u(t)\\ 2 < SUp \\u E (t) - Wool^ + SUp \\u(t) - Uool^ 
t>T t>T t>T 

= \\U £ (T) - MooHa + SUp \\u(t) - Moo^ 

t>T 

< \\u £ (T) - u(T)\\ 2 + 2 sup \\u(t) - M00II2, 

t>T 

hence 

sup \\u £ (t) - u(t)\\ 2 = 
t>o 

< 

Letting e — > + , the first term tends to because of the convergence result in [0,T]. 
Letting T — > +00, the second term tends to because u(t) — > u^. This proves (2.8) 
provided that initial data satisfy (4.3). 

Convergence for general data Let now uq £ — > uq be any family satisfying (2.7). Let 
us choose a sequence {won} Q L°°(Q) n W 1,00 (Q) with -u „ — > uq. For every n G N, let 
u e,n(t) be the solution of the approximating problem with «£,„(()) = u 0n , and let u n {t) 
be the solution of the limit problem with u n (0) = u 0n . For every nGNwe already know 
that u £)Tl — > u n in C° ([0, +00); L 2 (Q)), because in this case the sequence of initial data 
does not depend on e and satisfies (4.3). 

Since both the approximating problems and the limit problem generate a contraction 
semigroup in L 2 (Q), we have that 

\\u £ (t) - u(t) || 2 < \\u £ (t) - u e>n (t) || 2 + \\u £ , n (t) - u n (t) || 2 + \\u n {t) - u(t) || 2 

< ||w 0e - U 0n \\ 2 + \\u £ , n (t) - U n (t)\\ 2 + \\u 0n - U \\ 2 . 

Taking the supremum over all t > 0, and letting e — > + , we obtain that 
limsup sup \\u £ (t) - u(t)\\ 2 < 2\\uQ n - uq\\ 2 . 

e^0+ t>0 

Letting n — > +00, we finally obtain (2.8) for general data. □ 



max< sup \\u £ (t) — u(t)\\ 2 , sup \\u £ (t) — u(t)\\ 2 

Me[0,T] t>T 

sup \\u £ (t) - u(t)\\ 2 + 2 sup \\u(t) - Woolb- 

te[o,T] t>T 
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